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Abstract. Since gravitational wave spacetimes are time- varying vacuum solutions of 
Einstein's field equations, there is no unambiguous means to define their energy con- 
tent. However, Weber and Wheeler had demonstrated that they do impart energy to 
test particles. There have been various proposals to define the energy content but they 
have not met with great success. Here we propose a definition using "slightly broken" 
Noether symmetries. We check whether this definition is physically acceptable. The pro- 
cedure adopted is to appeal to "approximate symmetries" as defined in Lie analysis and 
use them in the limit of the exact symmetry holding. A problem is noted with the use 
of the proposal for plane-fronted gravitational waves. To attain a better understanding 
of the implications of this proposal we also use an artificially constructed time-varying 
non-vacuum metric and evaluate its Weyl and stress-energy tensors so as to obtain the 
gravitational and matter components separately and compare them with the energy con- 
tent obtained by our proposal. The procedure is also used for cylindrical gravitational 
wave solutions. The usefulness of the definition is demonstrated by the fact that it leads 
to a result on whether gravitational waves suffer self-damping. 
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1. Introduction 



Gravitational wave spacetimes are non-static vacuum solutions of the Einstein Field Equa- 
tions (EFEs), i.e. they have no timelike Killing Vectors (KVs) [IJ. This creates a problem 
with the definition of energy for gravitational waves in General Relativity (GR), as energy 
conservation is guaranteed for spacetimes that admit timelike KVs. Since for these waves 
the stress-energy tensor is zero, there was a debate whether they really exist [21 13]. To 
demonstrate their reality, Weber and Wheeler obtained (first and second order) approx- 
imate formulae for the momentum imparted to test particles in the path of cylindrical 
gravitational waves [H H] . Later Ehlers and Kundt did the same for plane gravitational 
waves [3]. Using the pseudo-Newtonian formalism Qadir and Sharif [5j obtained a gen- 
eral closed form expression for the momentum imparted to test particles in an arbitrary 
spacetime, which gave the Weber- Wheeler approximation for cylindrical waves. 

To define the energy content of gravitational waves, different people have attempted dif- 
ferent approximate symmetry approaches. One such attempt was to assume that conser- 
vation of energy holds asymptotically and examine whether it would work for gravitational 
waves assuming a positive definite energy [6] . An altogether different idea was adopted to 
provide a measure of the extent of break-down of symmetry by the integral of the square 
of the symmetrized derivative of a vector field divided by its mean square norm [71 [8]. 
This led to what was called an almost symmetric space and the corresponding vector field 
an almost Killing vector [9J. This measure of "non-symmetry" in a given direction was 
applied to the Taub cosmological solution [TD] and to study gravitational radiation. It 
provides a choice of gauge that makes calculations simpler and was used for this purpose 
Essentially based on the almost symmetry, the concept of an "approximate sym- 
metry group" was presented [12]. However none of them seem unequivocally successful. 
The approach of a "slightly broken symmetry" seems promising but merely providing 
simplicity of calculations is not physically convincing. Other approaches need to be tried 
to find one that is significantly better than the others, in that it is consistent with other 
physical concepts and leads to new physical insights. It was speculated that the use of 
approximate Lie symmetry methods for differential equations (DEs) [131 E] ^lay give 
a resolution to the problem of energy in non-static spacetimes [T^l US]- In this paper 
we apply these methods to propose a resolution of the problem of the energy content of 
gravitational wave spacetimes. 

By virtue of Noether's theorem [T7| for every infinitesimal generator of symmetry of a 
Lagrangian (called a Noether symmetry), there is a conserved quantity. For time trans- 
lational invariance it is energy that is conserved. It is for this reason that it was hoped 
that the symmetry approach could prove fruitful for defining the energy content of gravi- 
tational waves. Reducing from the maximally symmetric Minkowski spacetime (10 KVs) 
to the Schwarzschild [18] and Reissner-Nordstrom (RN) spacetimes linear and spin 
angular momentum conservation are lost. Using approximate Lie symmetry methods for 
DEs these conservation laws were recovered as trivial first-order and second-order ap- 
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proximate conservation laws respectively. Reducing from the Schwarzschild spacetime to 
the Kerr (or charged-Kerr) spacetime we lose angular momentum conservation. Going 
directly from the Minkowski spacetime to the Kerr (or charged-Kerr) spacetime we lose 
linear and spin angular momentum conservation as well. These lost conservation laws were 
recovered as "trivial" first-order and second-order approximate conservation laws for this 
spacetime [TB] (in that there is no exact symmetry part mixed in with the approximate 
symmetry infinitesimal generator). 

The trivial first order approximate symmetries did not provide any new insights. However, 
for the second-order approximate symmetries of the geodesic equations for the RN |T5] 
and charged-Kerr spacetimes the time translational approximate symmetry generator 
had to pick up a rescaling factor to provide a necessary cancelation. This corresponds 
to rescaling the energy of test particles and hence could possibly lead to a definition of 
energy in the spacetime. Clearly, the energy content would then be defined as the scaling 
factor. 

To check the proposal for defining energy in gravitational wave spacetimes by using slightly 
broken symmetry as defined in Lie analysis [13j, plane-fronted parallel-rays (pp) gravita- 
tional waves [19] are first investigated. For the approximate symmetries of pp-waves first 
the t-dependent part of this spacetime is removed to make it static and taken as the 
unperturbed spacetime. Then the exact pp-wave is taken as a perturbation on this static 
spacetime by considering the arbitrary amplitude of the wave as a small parameter, e. 
Since does not appear in the geodesic equations for perturbed pp-waves, there is a 
problem in applying the definition of second-order approximate symmetries of ordinary 
differential equations (ODEs), which gives the scaling factor mentioned earlier, to them. 
It can also be seen from the geometry of pp-waves in which the wave fronts are like moving 
parallel planes and the curvature of the spacetime is absolutely zero before the pp-wave 
pulse arrives and after it has passed pQ. There is no region where there is a slight shift 
from the flat geometry as required for obtaining an approximate symmetry. Thus the 
proposal for determining the energy content of pp-waves cannot be checked. The con- 
formally invariant Weyl tensor |20] represents a pure gravitational field. In some sense 
it tells us about the gravitational energy of the spacetime, but it does not give a direct 
measure of the gravitational energy. The stress-energy tensor gives the matter content 
of the spacetime [T]. For this perturbed spacetime the stress-energy tensor is zero while 
the Weyl tensor is nonzero. To obtain a better understanding of the energy rescaling 
in plane gravitational waves, the artificially constructed example of a plane symmetric 
"wave-like" spacetime [21], which represents a gravitational wave interacting with matter 
is investigated. Here we do obtain a scaling factor which gives the rescaling of energy in 
the spacetime field. 

For the plane wave-like spacetime [2T], along with trivial approximate symmetries, a non- 
trivial first-order approximate symmetry was found. The first-order approximate (stable) 
first integral corresponding to the non-trivial first-order approximate symmetry of this 
wave-like spacetime is calculated here. The first-order non-trivial approximate time-like 
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Noether symmetry is used with the momentum vector to obtain a conserved quantity 
which gives the energy non-conservation due to time variation. To check this quantity 
first-order approximate Noether symmetries of pp-wave spacetime are investigated. Only 
the exact symmetries are recovered as trivial first-order approximate Noether symmetries, 
which gives the exact conservation laws as trivial first-order approximate conservation 
laws. 

Next cylindrically symmetric exact gravitational waves [22], which are physically easier 
to understand, are considered. To study their approximate symmetries the t-dependent 
part is first removed to define a static spacetime and the exact wave is dealt with as a 
perturbation of this static spacetime, taking the strength of the wave as a small parame- 
ter, e. A scaling factor is obtained for this spacetime which gives the rescaling of energy. 
Noether symmetries for the cylindrically symmetric case are also considered to look at the 
conserved quantities. First a cylindrically symmetric wave-like spacetime is investigated 
which has a non-trivial first-order approximate Noether symmetry that gives the conserved 
quantity like the plane symmetric case. There is no non-trivial approximate symmetry 
for the perturbed cylindrical wave spacetime. The approximate Weyl and stress-energy 
tensors (up to first-order in e) for the cylindrical wave spacetimes are nonzero. Electro- 
magnetic waves in their interaction with matter get damped. This is known as Landau 
damping |23] . Since GR is a non-linear theory, gravitational waves have self-interactions. 
The question arises whether we should expect Landau self-damping. With our proposal, 
a self-damping is seen for cylindrical waves. 

The plan of the paper is as follows. The next section briefly reviews the mathematical 
formalism to be used. In section 3, second-order approximate symmetries of the geodesic 
equations for the plane wave spacetimes are discussed and graphs of the scaling factor for 
them are given. In the same section the Weyl and stress-energy tensors for plane wave 
spacetimes are also discussed. Noether symmetries of pp-waves are studied and a review 
of the approximate Noether symmetries for the plane symmetric wave-like spacetime is 
given, in section 4. In the next section, second-order approximate symmetries of the 
geodesic equations for the cylindrically symmetric case are provided and the graphs of 
the scaling factor for them given. The Weyl and stress-energy tensors for the cylindrical 
waves are also discussed in the same section. In section 6 approximate Noether symmetries 
for the cylindrical wave spacetimes are investigated. Finally a summary and discussion 
are presented in section 7. 



2. Review of mathematical formalism used 

We first review the definition of the second-order approximate symmetries of a system of 
ODEs under point symmetries. A vector field 




(1) 
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is called a second-order approximate symmetry of the system of perturbed ODEs 

E = Eo + eEi + e'E2 + 0(e=^), (2) 
if the following condition ([12] and references given there in) holds 

(Xo + eXi + e^X,) (Eo + eE, + e'M^^^^^^^^,.^^^^^,s^ = 0{e% (3) 
where Xq is the exact symmetry generator of the system of ODEs Eq, i.e. 

(Xo)(Eo)|e„=o = 0, (4) 

Xi, X2 are the first-order and second-order approximate parts of the approximate sym- 
metry generator respectively, Ei is the first-order perturbed part and E2 is the second 
order perturbed part of the system of ODEs respectively. It should be noted that the 
scaling factor comes from the applications of the perturbed system of DEs in subscript of 
(3), as required. 

Noether symmetries are those infinitesimal symmetry generators that leave a Lagrangian 
L{s,x^,x^) invariant. They form a Lie algebra that contains the isometrics for the La- 
grangian that minimizes arc length, with at least one extra symmetry, d/ds, [21]. It is 
defined vector field [H EHl ESI 



X=e(^>a:-)^ + ^^(^,a;-)^, (5) 
where /x, z/ = 0, 1, 2, 3, such that 

xWl + (Z},o^ = ^.A (6) 

where 74(s,x'^), is a gauge function. The total derivative operator and the first pro- 
longation Xt^l of the vector field X given by (5) are 

and 

xw = X + (r^;; + - ^,sx^ - 

For more general considerations see [14j . The significance of Noether symmetries is clear 
from the following theorem [T7], proved in [27] . 

Theorem 1. If X is a Noether point symmetry corresponding to a Lagrangian L{s, x^, x^) 
of a second-order ODE x^ = g{s, x, x^), then 

I = iL + {r^^-x^i)—-A, (9) 



is a first integral of tlie ODE associated witli X. 



First-order approximate symmetries of a Lagrangian (or first-order approximate Noether 
symmetries) [261 [2H] are defined as follows. For a first-order perturbed system of ODEs 

E = Eo + eEi =0(6^), (10) 

corresponding to a first-order Lagrangian, which is perturbed up to first-order in e, 

L(s, xf", x^, e) = Lo(s, x^ , x^) + eLi(s, x^ , x^) + 0(6^), (11) 

the functional JyLds is invariant under the one-parameter group of transformations with 
approximate Lie symmetry generator 

X = Xo + eXi + 0(e2), (12) 

up to gauge 

A = Ao + eAi, (13) 

where 
if 

X^Lo + {Ds^o)Lo = DsAo, (15) 

and 

X^Lo + X^Li + (D,ei)^o + {Ds^o)Li = D.,A^. (16) 

Here Lq is the exact Lagrangian corresponding to the exact equations and Lq + eLi the 
first-order perturbed Lagrangian corresponding to the first-order perturbed equations. 
The perturbed equations (3) and (16) always have the approximate symmetry generators 
eXo which are known as trivial approximate symmetries and X given by (1) and (12) with 
Xq 7^ is called a non-trivial approximate symmetry. These approximate symmetries of 
a manifold, form an approximate Lie algebra |29j . 

The first-order approximate first integrals are defined by setting I = lo + eli, L = Lo + eLi, 
^ = Co + ^^1, V = Vo + ^Vi and A = Aq + eAi in (9) and equating the coefficients of like 
powers of e on both sides. This gives the zeroth (exact part) and first-order approximate 
part of the first-order approximate first integrals 

/o = ^0^0 + ivH - ^%)^ - ^0, (17) 

h = CoL^ + C^Lo + H - i-eo)^ + « - i'^G)^ - A,. (18) 

If Jo vanishes, then / is called an unstable approximate first integral and is otherwise 
called stable. A detailed discussion on the approximate first integrals for Hamiltonian 
dynamical systems is given in [30] . 



6 



The Weyl tensor in component form is given by 



C\cd = R'^hcd - 2^^cRbd - S^Rbc + QbdR^c - QbcR^d) + QR{^d9bc - S'^gbd)- (19) 

Here R^'fy^d is the Riemann curvature tensor, Rat is the Ricci tensor, R is the Ricci scalar, 
Qab is the metric tensor and 6^ is the Kroneker delta. For a 4 dimensional spacetime 
the Weyl tensor has 10 independent components [31]. If the Weyl tensor vanishes in a 
neighborhood of a spacetime, the neighborhood is locally conformally equivalent to the 
Minkowski spacetime. Thus the Weyl tensor has geometric meaning independent of any 
physical interpretation. 

The stress-energy tensor can be calculated from the EFEs 

Tab=-{Rab-lRgab), (20) 

where k is the gravitational coupling. For a 4-dimensional spacetime this tensor has 
10 independent components. At each event of the spacetime this tensor gives the energy 
density, momentum density and stress as measured by any and all observers at that event. 
Since for gravitational wave spacetimes Tab is always zero and C°^^^ is nonzero, there is 
no matter or energy or momentum however there is the Weyl curvature. If there is no 
mass or energy at a given event, the Ricci tensor vanishes. If it were not for the Weyl 
tensor, this would mean that matter at one place could not have gravitational influence 
on distant matter separated by a void. Thus the Weyl tensor represents that part of 
spacetime curvature which can propagate across and curve up a void. 



3. Second-order approximate symmetries and energy 
rescaling: plane wave spacetimes 

The line element for pp- waves [12], is 

ds^ = huj^[{x^ - 2/^) sm{uj{t - z)) + 2xycos{uj{t - z))]{dt^ + dz^ - 2dtdz) + dt^ 

~dx^ - dy^ - dz^, (21) 

where h is the amplitude of the wave and u is the frequency. 

Now we remove the t-dependent part of (21) and putting /i = 1 to define a static spacetime 

ds"^ = uj^lix^ - y2) + 2xy]{dt'^ + dz'^ - 2dtdz) + dt^ - dx^ - dy^ - dz'^. (22) 

To obtain the approximate symmetries of pp-waves the exact pp-waves are considered as 
a perturbation on the static spacetime (22). For this purpose the amplitude /i = e, is 
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taken as a small parameter and the line element of the perturbed pp-waves is 



t) 




t) 




0, 


(24) 


0, 


(25) 


0, 


(26) 


t) 




t) 




0. 


(27) 



ds'^ = uj'^[x^ -y^ + 2xy + e{(x^ - y^) sin(cj(t - z)) + 2xycos{uj{t - z))}]{dt'^ 

+dz^ - 2dtdz) + dt^ - dx^ - dy^ - dz^. (23) 

For this perturbed pp-wave spacetime (23) we have the system of first-order perturbed 
geodesic equations and there do not appear e^, 

^3 

t + uj'^it — z){{x + y)x + (x — y)y} + e[— {(x^ — y'^) coiiuj{z 

+2xysm.uj{z — t)}{p + — iz) — u'^it — z){x sinco'(z 
—ycosuj{z — t)}x + uj'^{ysm.uj{z — t) + x cosuj{z — t)}y\ ■ 
X + [uj'^{x + y) — euj'^{xsmuj{z — t) — ycosuj{z — t)}](P + — iz) ■ 
y + [uj'^{x + y) — euj'^{x cos uj{z — t) — ysmiu{z — t)}]{P + z^ — iz) ■ 

z + uj'^{i — z){{x + y)x + (x — y)y} + e[— {(x^ — y'^) cosuj{z 

-|-2x?/sinci;(2; — t)}{P + — iz) — uj'^{i — z){x sm.uj{z 
—ycosuj{z — t)}x + uj'^{y smuj{z — t) + x cosuj{z — t)}y] ■ 

Since there is no quadratic term in e, in the above geodesic equations, we cannot apply the 
definition of second-order approximate symmetries, which gives us the energy rescaling 
factor to them. This behavior is consistent with the pp-wave geometry in which the wave 
front moves as parallel planes and the spacetime curvature is zero before and after the 
pp-wave pulse pQ. 

To obtain a better understanding of the energy rescaling in plane gravitational waves we 
apply the definition of second-order approximate symmetries of ODEs to the second-order 
perturbed geodesic equations (29) - (32) for the artificially constructed example of plane 
symmetric wave-like spacetime [21] • For this purpose a non-fiat plane symmetric static 
spacetime [32], was considered with /i(x) = z/^(x) = (x/X)^, 

ds^ = e^'^'-^^dt^ - dx^ - e2'^(")(di/2 + dz^), (28) 

where X is a constant having the same dimensions as x. Since gravitational waves are 
non-static spacetimes therefore the static spacetime (28) was perturbed with a time- 
dependent small parameter (for definiteness by et) to make it slightly non-static. For this 
the metric (28) was taken with z/(x) = {x/X + ^t/T) and /i(x) = (x^/X^ -|- et/T), where 
T is a constant having dimensions of t. Retaining and neglecting its higher powers 
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second-order perturbed geodesic equations are obtained 

i+ |tx - ^[t^ - if + i2)e2((^/^)^-/^)] + ^[P + {f+ 

i^y(i-/x?-m^+O{e^) = 0, (29) 

+i2)g2(./x)^] + |^[i'e^^'/^ - ^{f + i2)e2W^)^] + 0{e') = 0, (30) 

A i ^ X 

Ax 2e 2te^ 
i/ + —xy + -ty - —ty + 0{e') = 0, (31) 

Ax 2e . 2te^ ■ 
z + —xz + -tz - —tz + 0{e^) = 0. (32) 

The Lie symmetry algebra of the exact or unperturbed geodesic equations (i.e. when e = 0, 
in (29) - (32)) includes the generators of the dilation algebra d/ds, sd/ds, corresponding 
to 

C{s) = CqS + ci. (33) 

In the determining equations for the first-order approximate symmetries of the geodesic 
equations for the Schwarzschild spacetime [18] the terms involving C,s = cq cancel out. 
Here the terms involving C,s = Cq do not cancel automatically but, like the RN [15] and 
charged-Kerr [TH] spacetimes, collect a scaling factor to cancel out. In this case the scaling 
factor is 

1 ^ 

Energy conservation is related with time translation and ^ is the coefficient of djds in 
the point transformation given by (5), where s is the proper time. The scaling factor 
(34) corresponds to the rescaling of energy of a test particle in this wave-like spacetime 
field. Since the scaling factor (34) involve the derivatives of the coordinates and the 
derivatives only apply to the paths of the particles. To get energy in the spacetime field 
the derivatives of the coordinates are replaced by the first integrals. Therefore we get 

^ [e-^^/^ + 2e-2(-/^)(-/^+i)] . (35) 

This energy expression is plotted below for different values of t and x, using Mathematica 
5.0. The values of X and T are arbitrary. The above scaling factor for this wave-like 
spacetime depends linearly on t and in both diagrams below the energy in the gravitational 
field increases linearly with time. In Fig. 1 the energy is seen to decrease along x and 
disappear sharply close to a; = 0. To see the variation with x we enlarge the diagram by 
reducing the range of x in Fig. 2. As we move along x the increase in energy with time 
becomes gradual. Since the small parameter e, (which is considered as the strength of 
the wave) is arbitrary the units of energy are arbitrarily chosen. Throughout this paper 
gravitational units are used and space, time and mass are given in seconds. 
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Figure 1: Plane symmetric gravitational wave-like spacetime. The energy increases in- 
definitely in time close to x = and then disappears suddenly after some distance. The 
small parameter e, (considered as strength of the wave) is arbitrary in all the spacetimes 
discussed in this paper. Thus the units of energy are chosen arbitrarily. Throughout this 
paper gravitational units are adopted and space, time and mass are given in seconds. 




Figure 2: This is an expanded version of Fig. 1. Here the range of x is shrunk and it is 
seen that the energy decreases smoothly with distance. 
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Though the Weyl tensor gives information about the gravitational energy of the spacetime, 
it is not clear how to obtain a measure of the energy in it. For the pure gravitational part 
of the perturbed pp-wave spacetime the independent nonzero components of the Weyl 
tensor are 

-uj'^[uj'^{x'^ — y'^ — 2xy) + 1 + e{2u'^(x'^ — v"^ — xy) sinu;(2; — t) — 2uj'^xy 
t) - smu;{z - t)}] + 0{e^) = C\,, = C\,, = -C\, = C\,, = C\,,, 
-uj'^[uj^{y'^ -x^ - 2xy) + 1 + t{uj'^{x^ - y"^) sincj(2 - t) - - 

Axy) cosuj{z -t) + coscu(2 - t)}] + 0{t^) = C\,, = = C\,,. (36) 



'-^ 101 - 
COSCo'(z 

102 " 



In the literature [33] the Weyl tensor is usually defined with valence (1, 3). In spinors it 
is naturally given as a tensor of valence (0, 4) [M]. For usual purposes the form does not 
matter, but for differential symmetries of the tensor the form is crucial [35]. In covariant 
form the components of the Weyl tensor are 

Coioi = — — esinci;(2; — t)] = C0113 = C1313 = — C0202 = C0223 = (^2323; 

C0102 = -uj^[l + esmu{z - t)] = C0123 = C0213 = C'i323- (37) 

From here it appears that the (0, 4) form may give the physically relevant quantities as 
the space dependence in (36) does not seem to correspond to the geometry of the pp- 
wave, while (37) does. Here the pure gravitational field which "curves up the void" is 
seem to be sinusoidal. For this spacetime there is obviously no nonzero component of the 
stress-energy tensor. 

There are 6 nonzero components (up to first-order in e) of the Weyl tensor for the above 
plane symmetric wave-like spacetime 

C%i = ^(2x + X) + 0(e^), 

C\o2 = = -^(1 + e^)(2x + X) + O(e^), 

C 212 — ^ 313 ~ ~^°2025 ^ 323 ~ ~2C%2- (38) 

From Figs. 1 and 2 (where the wave is along the x direction), it is clear that the energy 
in the gravitational field of the plane wave-like spacetime increases with time. Therefore 
the first component of the Weyl tensor must depend on t linearly which corresponds to 
the covariant form (given below) and not the mixed form. 

2x/X Of 

C^oioi = + e^)(2^ + + 0{e% 

2x'2/X'2+2x/X Af 

C,2m = Co303 = (1 + e-)(2x + X) + 0{e'), 

C1212 = C1313 = —(^0202, C2323 = — 2Co202- (39) 
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The nonzero components of the stress-energy tensor for this wave-hke spacetime are 
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T22 = T33 = -^^a + e^)(2xX + Ax' + 3X') + 0(e2), 

To, = ^,{x + X) + Oie'). (40) 

It is worth noting that the x-direction stress has no approximate part of first-order and 
the approximate part of the energy increases hnearly with time and quadratically at large 
distances. More interestingly there is an approximate momentum in the x-direction that 
increases linearly with the value of x. This linear increase in energy was built into the 
metric and it entails the momentum in the x-direction. 

We give the ratio of energy density imparted to the matter field 

where {Tqq)e and (Too)p are the energy densities of the exact (i.e. when e = 0) and first- 
order approximate spacetimes respectively. For the plane symmetric wave-like spacetime 
we have 



4. Approximate Noether symmetries of plane wave 
spacetimes 

The Lagrangian defined for (28) is |2T] (throughout this paper the Lagrangian of a space- 
time would means the Lagrangian for the geodesic equations of the spacetime) 

L = e'^/^'P -x'- e'^'/'^'iy' + z'). (43) 
Its symmetry generators are 

d d d d d d 

Xo = ^, Xi = — , X2 = — , -K-i = y—-z — , Yo = — , A = c, (44) 
at oy oz oz ay as 

where c is a constant, Xq corresponds to energy conservation, Xi and X2 correspond to 
linear momentum conservation along y and z, while X3 corresponds to angular momentum 
conservation in the yz plane [36i| . 
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The first-order perturbed Lagrangian is 

L = e^^/^i^ - - e^^'/^'(y^ + i^) + ^^[e^^/^i^ - e^^'/^'(y^ + i^)] + O(e^), (45) 

yielding the non-trivial approximate symmetry 

d \ d d d 

along with the trivial symmetries, and the gauge function Ai is again a constant. The 
stable first integral for the symmetry generator (46) is 

/ ^ 2e2-/^i + ^ [e^^/^tt + e2(-/^)' {yy + zz)]. (47) 

Contract the energy momentum vector with the time-like approximate Noether symmetry 
generator (46), to obtain the conserved quantity 

Q^E-^{tE + ypy + zp,), (48) 

where E is the energy and p is the momentum. This gives the energy non-conservation 
due to time variation. That is the energy imparted to the test particles with energy and 
momentum given by (48). However this does not give the energy in the spacetime field. 

To check the conserved quantities in the pp-wave spacetime, we investigate the first-order 
approximate Noether symmetries for this spacetime. The Lagrangian for the pp-wave 
spacetime (21) is 

L = huj^[{x^ - y^) sm{u{t - z)) + 2xy cos(a;(t - z))] (P + - 2iz) +f-x'^-f-z^. (49) 
This Lagrangian admits the following symmetry generators 

d d d 
Xo =7^ + 7^, Yo = — and A = c (constant). (50) 
at az as 



The Lagrangian for the static spacetime (22) 

L = uj^[{x^ - 7/2) + 2xy]{i^ + - 2iz) + f - x^ - f - (51) 
has 3 symmetry generators 

d d d 

^'^d~r ^'^d~z' ^'^ds' ^^^^ 

and the gauge function is a constant. 
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The Lagrangian for the perturbed pp-wave spacetime (23) is 

L = u?{x^ - 2/^ + 2xy + e{{x^ - y^) sin(u;(t - z)) + 2xyco^{yj(t - z)))\(p + 

-2iz) + - - ij"^ - z^ . (53) 

For e = 0, the above Lagrangian (53) reduces to (51). Using this first-order perturbed 
Lagrangian and the three exact symmetry generators given by (52) in (16), in the resulting 
system of determining equations two constants corresponding to the exact symmetry 
generators appears. These two generators have to be ehminated for consistency of the 
determining equations, making them homogeneous. The resulting system is the same as 
for the static (exact) spacetime (23). Thus there is no non-trivial approximate symmetry 
for this perturbed Lagrangian and the gauge function is a constant. Hence we can not 
obtain the conserved quantity in the case of perturbed pp-wave spacetime. Only the 
three exact symmetry generators are recovered as trivial first-order approximate Noether 
symmetries which gives trivial first-order approximate conservation laws for energy and 
linear momentum along z. 



5. Second-order approximate symmetries and energy 
rescaling: cylindrical wave spacetimes 

The line element of the cylindrically symmetric exact waves [22] is 

ds^ = e'^''-^\df - dp') - p'e-'^dct)'' - e'^dz\ (54) 
where 7 and if) are arbitrary functions of t and p, subject to the vacuum EFEs 

^"+V-V^ = 0, y = p(7/;'2 + V''), ^ = 2p^^\ (55) 
P 

where dot denotes differentiation with respect to t and prime with respect to p. The 
solution of (55) is given by [1] 

ip = AJq{up) cos{ujt) + BYq{up) sin(ct;t), (56) 

7 = ^Ljp[{A'JoJo' - B'YoYo') cos(2cut) - AB{{JoYo' + YqJo') sm{2ut) 

-2{JoYo' -YoJo')ujt}]. (57) 

This metric has two KVs d/dcp and d/dz [3ij; this means that there is only azimuthal 
angular momentum conservation and linear momentum conservation along z. 

To discuss the approximate symmetries of cylindrical waves first a static spacetime is 
defined as follows. We remove the t-dependent part in (54) and put the strength of the 
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wave A = 1 and B = 0. Yq become singular at the origin and putting B = give us rid 
of this singular behavior. 

where 

ipo = Joi^jp), 7o = —Jo{ujp)J'q{ujp). (59) 

For the approximate case we put the strength of the wave as a small parameter i.e. A = e 
and take the exact wave as a perturbation on the static metric (58) in the following way. 

^p = Jo{ujp){l + ecos{ujt)) = ipo + e^pi 

7 = '^Mup)J'o{up){l + cos{2ut)) = 70 + e^- (60) 

Thus the second-order perturbed geodesic equations are 

t + 2(7^ - ij'o)ip - e[Mt^ + p2) + p2e-27o^^02 _ e2{2Vo-7o)^^^-2 + 2ij[ip] + e^lii 

(t2 + p2) + 4e2(2^o-^oViV^i-2'^ - 27; tp] + O(e^) = 0, (61) 
P+il'o- i'oW + + ip) + pe-2^»(^[, - 1)02 + e2(2V'o-7o)^/^2 _ ^[^/ ^ ^2) 
-p2e-2^°V^i02 - e2(2V^«-^«)(4V^^,V^i - V^'Ji^ + ^;ip] + e^Yi[{t^ + p2) _ p^^e-^i^{^'^ 

-1)02 + 2e2(2'^°-^°)(4V';V'i' - li-^'o - + l[ip] + 0(6^) = 0, (62) 

4>+kl- ^P',)p<P - e{i;[p + - ^(1 - ^o)p0 + Oie') = 0, (63) 

z + ^'opz + e{i)[p + i)ii)z - 2eV[,V'2pi + O(e^) = 0. (64) 

The dot and prime over 70, ipo, 71 and ipi denote derivatives with respect to ojt and ojp 
respectively. For this perturbed wave spacetime the scaling factor is 

iiit^ + p2) + 4e2(2V'o-7o)^^^^^-2 _ 2^[ip. (65) 

To replace the derivative of the coordinates t, z and p the exact first integrals and the 
metric (58) are used. Further it is assumed that there is no initial velocity in the z and 
directions. Hence i and vanishes and the following scaling factor is obtained 

^_^g2(i/'o-7o) j'g2{?/.o-7o) _|_ g3(i/)o-7o) — 1] — 2'y[e^^''^°~'^'^\e^^''^°~'^'^^ — l-Y^"^, (66) 

where 71 is given in (60). This scaling factor involves the Bessel function of first kind 
and its derivatives. The asymptotic representation of the Bessel function of first kind for 
large value of the argument is given in [SS]. Using this asymptotic representation of the 
Bessel function in (66), we obtain an asymptotic representation of it as follows 

^^|^[(|cos(.;p)|)3/2sin(2.;t)](^p)-i/2 + 0(M-3/2). (57) 

Thus the energy in this perturbed spacetime field is rescaled by the factor (67). It is 
plotted below for different values of t, p and uj (in radians per second), in which the 
energy oscillates between positive and negative values and goes to zero as p tends to 
infinity. Here the behavior is much more recognizably wave-like. Since the strength of the 
wave, A = e, is arbitrary the energy is given in arbitrarily chosen units. 
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Figure 3: Cylindrically symmetric gravitational waves with uj = 15. The gravitational 
energy oscillates between positive and negative values and disappears as p approaches very 
large value. The units of energy are arbitrary in all diagrams. 




Figure 4: To see the behavior of energy for comparatively larger distance, therefore the 
range of p is extended to 100 units. 
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Figure 5: To see a further extended version of the above fig. 4; the rang of p is given in 
units 0/ 10^. 




Figure 6: Here the value of the frequency is comparatively small i.e. 00 = 0.05. To see the 
variation along time, therefore the rang oft is kept larger. 
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For completeness we also investigate second-order approximate symmetries of the geodesic 
equations for the cylindrically symmetric wave-like spacetime. For this purpose a cylin- 
drically symmetric static metric is taken 



ds^ = e'^p'^dt^ - dp" - e^''^p\a'd(l? + dz"), (68) 

with v{p) = {p/ RY and p{p) = (p/R)^, where i? is a constant having the same dimensions 
as p. For the approximate symmetries of this cylindrical wave-like spacetime, z/(a;) = 
{p/Ry + 2et/T and p{x) = (p/R)^ + 2et/T are taken in the metric (68), where T is a 
constant having dimensions of t. Like the plane symmetric case the scaling factor for this 
cylindrical wave-like spacetime is 



4T2 



Thus the energy in the field of this wave-like spacetimes is rescaled by the factor (69). The 
plots for this case, are similar to those for the plane wave-like case discussed in section 3, 
with X replaced by p. 

The nonzero components of the Weyl tensor for the cylindrically symmetric perturbed 
wave spacetime are 

C\, = - 7o + 2^;' - + <^fo€ + - i^i - 2^; + 4^;^i 

-4V;(,Vi - 2<V^i + 27;Vi)] + 0(e2), 

C\2 = jl^i^o - 7o + 8^o' + H - 2^; - Qi^Wo + + ^1 + SV^oV^'i 

C°303 = ^[2^0 + 7o + 4^o' + 37o + 2^o - Q^Wo + H^i + 2^i 

+47A;^; + ^;-37;^'i)] + 0(e2), 

^1 _ _ ~4ijQ^0 ^1 _ i 2{2Vo-7o)(^0 

212 ~ 303' ^ 313 ~ ^ 202 

2(2,/.0-70) .. 1 

C\,, = [7^' - < + ij'o- 2^;^ - 2e{fl - ^1 + 2^;^; - -^0] + 0(6^), 

6 p 

C\, = -6p2e-2^°[(2^; + 7;)^i + A) + 0(6^), 

C°3i3 = 6e^(^^°-^«)[(2^; - 7;)V^i + ^i] + 0{e'). (70) 

This yields the pure gravitational field for this cylindrically symmetric perturbed wave 
spacetime. As it is evident from Figs. 3 to 6, that the energy in the gravitational field 
oscillates and then vanishes for large p, here all the components of the Weyl tensor also 
depend on the Bessel function of the first kind and its derivatives which oscillates and 
goes to zero as p approaches very large value. The last two components only appear for 
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the approximate part of the spacetime. In this case the components of the Weyl tensor in 
the covariant form are not very different from those in the mixed form given above and 
therefore we do not give them separately. 

The non-vanishing components of the stress-energy tensor are 

Too = Tn = k[< - 7^ + 2ei^'M + ^(e^), 

(V'o - 7o - ^0' + SV'^, - i^'oi^', + 2(V^i - + 0(e2), 

Toi^^eV'oV^i + Ole^). (71) 

Like the components of the Weyl tensor the above components of the stress-energy tensor 
also depend on the Bcsscl function of the first-kind and its derivatives. In this case of 
cylindrical perturbed waves the fraction of energy density imparted to the matter field is 

W — Or "^o"^! 

-C/imp — ^fc ,,2 „ 7 • \'^) 
YO 10 



The non- vanishing components of the Weyl tensor for the cylindrical wave-like spacetime 
are 

C\oi = ^(3^V - 2 V + - R') + O(e^), 

C%2 = a'C%,, = -(1 + 6|)^^(3i?V - 2V + 3/ - R') + 0{e% 

C 212 — ^ 313 ~ ~^°2025 ^ 323 ~ ~2C°202- (73) 

In covariant form the components of the Weyl tensor are 
<^oioi = ^(1 + e|)(3i?V - 2 V + 3p' - R') + 0(6^), 
C0202 = a^Coaos = -^^^"^^^ (1 + e^)i^R'p - 2Rp' + 3p' - R') + O(e^), 

C1212 = Ct^C'l313 = ~C'o2025 ^^2323 = ~2Co202- (74) 

The components in the covariant form arc physically reasonable as they follow the geom- 
etry of the constructed metric and the energy defined by approximate symmetry. 
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The nonzero components of stress-energy tensor are 

roi = e^(3p-2i?) + 0(62). (75) 

Here the momentum density is along radius of the cyhnder. For this case we have the 
same relative energy density imparted to the matter field, as given by (42). 



6. Approximate Noether symmetries of the 
cylindrical wave spacetimes 

The Lagrangian of the spacetime (68) is 

yields the following symmetry generators 

= = = = ^ = 

where c is a constant, Xq corresponds to energy conservation, Xi corresponds to azimuthal 
angular momentum conservation and X2 to linear momentum conservation along z, while 
X3 corresponds to angular momentum conservation. 

The first-order perturbed Lagrangian for the cylindrical wave-like spacetime is 

L = e^p/R?i^ - p2 _ eW«)^>2^2 ^ -2) ^ '^^^{p/R?i2 _ ^(p/R)\a^'' + z^)] + 0{^). (78) 

For this Lagrangian along with the exact symmetry generators given by (77), the non- 
trivial approximate symmetry Xo given by (79) is obtained. The gauge function Ai is 
again a constant, 

d \ d d d 

The corresponding stable first integral is 

I = 2e(''/^)'i + ^[e(^/^)'tt + e(^/^)'(a2# + zz)]. (80) 
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As for the plane wave-like spacetime the following conserved quantity corresponding to 
(79) is calculated 

Q^E-^{tE + (t>p^ + zp,) . (81) 



Using the Lagrangian for the spacetime (54) 

in (6) and solving the resulting system of determining equations we obtain the symmetry 
generator d/ds along with the two KVs d/dcj), d/dz and the gauge function A is a constant. 

The Lagrangian for the spacetime (58) is 

L = e^(70-Vo)(^2 _ ^2) _ ^2g-2V,o^2 _ ^2^,-2^ (gg) 

which admits the following 4 symmetry generators along with the gauge function as a 
constant 

d d d d 



The first-order perturbed Lagrangian for the cylindrical wave spacetime (54) with ip and 
7 defined by (60) is given by 

L ^ e^{jO-fo)^2 _ .2^ _ p2g-2V-o^2 _ g2V-o^2 _ 2e^i [c^Cto-V-o) (^2 _ ^2) _ ^2g-2V.o02 

+e''^°i']+0(e'). (85) 

For e = this Lagrangian reduces to the Lagrangian (83). Using this perturbed La- 
grangian and the 4 dimensional exact symmetry algebra of the static spacetime (58) in 
(16), we get a set of determining equations in which only 1 constant corresponding to the 
exact symmetry generator appears. This exact symmetry generator has to be eliminated 
for consistency of these determining equations, making them homogeneous. The result- 
ing system is once more the same as for the spacetime (58), yielding first-order trivial 
approximate symmetry generators. Thus there is no non-trivial approximate symmetry 
for this perturbed Lagrangian and the gauge function A is a constant. Hence energy con- 
servation, azimuthal angular momentum conservation and linear momentum conservation 
along the axis of the cylinder are obtained as trivial first-order approximate conservation 
laws. (Note that the technically "trivial" law may be physically non-trivial.) 



7. Summary and discussion 



The problem of energy in gravitational wave spacetimes using approximate Lie symmetry 
methods for DEs is addressed. To resolve this problem we used the second-order approx- 
imate symmetries of the geodesic equations for perturbed gravitational wave spacetimes 
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discussed here. First the pp-wave spacetime is investigated. Since there is no in the 
geodesic equations for the perturbed pp-waves, the definition of second-order approximate 
symmetries of ODEs which gives the scahng factor, can not be apphed to them. This is 
similar to the result of Qadir and Sharif's work [5], using the pseudo- Newtonian formal- 
ism, which just gave a constant momentum imparted to test particles in the path of the 
waves and no determinable value for it. For a better understanding of the implication 
of the definition of second-order approximate symmetries of ODEs, in plane symmetric 
waves this definition has applied to the artificially constructed time-varying non-vacuum 
plane symmetric spacetime [21], for which the scaling factor (35) is obtained. It is seen 
from the plots of the plane symmetric wave-like spacetime that the energy increases with 
time close to the origin for x and then disappears. Then we investigated the second-order 
approximate symmetries of the geodesic equations for the cylindrical wave spacetimes. 
The scaling factors (67) and (69) are obtained for these spacetimes. In the factor (67) the 
magnitude of the coefficient of (a;p)~^/^ is greater then the magnitude of the coefficient 
of (cup)"^/^, therefore, the contribution of the second term is very small and is neglected. 
This factor is plotted for different values of t, p and uj. It shows a behavior much more 
recognisably wave-like. In Figs. 3 to 6 the energy oscillates between positive and negative 
values along t and p. It disappears as p tends to infinity. 

To obtain the pure gravitational field and the matter field the approximate Weyl and 
stress-energy tensors for the gravitational wave spacetimes are calculated. The compo- 
nents of the Weyl tensor are given in the (0, 4) (covariant) form as well. For the perturbed 
pp-wave spacetime it appears that the (0, 4) form gives the physically relevant quantities 
as the space dependence in the (1, 3) (mixed) form of the Weyl tensor does not seem to 
correspond to the geometry of the pp-wave, while the covariant form does. For the wave- 
like spacetimes the components in the covariant form are physically reasonable as they 
follow the geometry of the constructed metrics and the energy defined by approximate 
symmetry. The stress-energy tensor density imparted to the matter field in the wave-like 
and perturbed cylindrical wave spacetimes was obtained. 

In GR, different people have tried to introduce the concept of a pseudo-tensor, to define 
energy and momentum. In this regard first Einstein obtained a pseudo-tensor to define 
energy in GR p9]. Following Einstein's idea, Landau-Lifshitz |10], Papapetrou [H] and 
Weinberg [52] gave different pseudo-tensors to represent the energy and momentum of the 
gravitational field. The idea, of introducing a pseudo-tensor has been criticized because 
all the pseudo-tensors are coordinate dependent and hence non-tens orial. This violates 
the basic principles of GR. Because of the coordinate dependence, many others, includ- 
ing MoUer [lU HI], Komar [15], Ashtekar-Hansen [^ and Penrose [17], have proposed 
coordinate independent definitions. MoUer realized that the use of a tetrad as the field 
variable, instead of a metric, makes it possible to introduce a first order scalar Lagrangian 
for the EFEs. Komar introduced a tensorial super-potential which is independent of any 
background structure and has uniqueness property. Ashtekar and Hansen defined the 
angular momentum in their specific conformal model of the spatial infinity as a certain 2- 
surface integral near infinity. Penrose defined quasi-local energy-momentum and angular 
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momentum using twistor-theoretical idea. However, each of these, has its own drawbacks 
im HHl So] . A detailed discussion on different definitions of gravitational energy is avail- 
able in "Quasi-local energy-momentum and angular momentum in GR: a review article" 
j50] and "Energy and momentum in GR" [51j- Using the idea of pseudo-tensors different 
people claimed that the gravitational energy should be positive at large scales as well as 
at small scales ESI ISH ESI ES] • The positivity of gravitational energy does not seems 
convincing because the total energy of the universe is zero [57], which suggests that the 
gravitational energy must fluctuate between positive and negative values to give the net 
energy of a spacetime zero. 

Our definition of gravitational energy, obtained from approximate Lie symmetry methods, 
avoids the pseudo-tensor and hence does not violate GR. Our expression of energy is 
also reasonable as the gravitational energy oscillates over positive and negative values, 
as it should. Admittedly, in the artificial example we constructed the energy increased 
linearly without limit. This was because of the (nonphysical) choice of a linearly increasing 
component of the metric tensor for convenience of computation, leading to a corresponding 
increase in the scaling factors (35) and (69). For the physical example of cylindrical exact 
waves, the Bessel function of the first kind goes to zero asymptotically for large values of 
the argument [38] . Correspondingly, our scaling factor (67) for cylindrical waves dies out 
asymptotically, giving a net zero energy. 

One of us (AQ) was drawn to study the energy of gravitational waves when a question 
was posed [SH] whether there is the analogue of Landau-damping of electromagnetic waves 
for gravitational waves. Since Maxwell's theory of electromagnetism is linear, electromag- 
netic waves do not interact with the field but are damped due to their interaction with 
matter. On the other hand GR is non-linear and so gravitational waves undergo self- 
interaction. This gives rise to the possibility of "Landau self-damping" of gravitational 
waves. On the other hand, the Khan-Penrose jSH] and Szekeres [20] solutions of colliding 
plane gravitational waves suggest that there could even be enhancement of the waves, as 
they lead to curvature singularities after the collision. The problem of definition of energy 
in GR makes it very difficult to answer the question posed. Using Wheeler's "poor man's 
approach", we can ask whether "the mass equivalent to the energy of the gravitational 
wave attracts and hence damps the waves" , or like the black hole, "the energy enhances 
the mass and hence the energy equivalent to it in the wave" . With our present proposal 
the question seems to be answerable. Classically the energy density in cylindrical waves 
reduces by the factor l/(2'7rp). From (67) the energy density decreases by a further factor 
of (3 X 2^'^^'^)/ ^y^^^^~^<llJf)j. Hence for sufficiently large p the scaling factor ~ 1/ ^/ujp^ 
which is a significant self- damping of the waves! This enhanced asymptotic attenuation 
of gravitational waves will obviously have profound observational significance. 

It would be of great interest to apply this approximate symmetry analysis to the Khan- 
Penrose and Szekeres solutions to see whether they suffer self-damping or enhancement 
according to our definition. Of course, it may be that the procedure will be inapplicable 
for those plane wave solutions as well. Also, the analysis should be applied to "spherical 
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solutions" like those of Nutku . 
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